arXiv:1506.08849v2 [hep-ex] 20 Aug 2015 
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Differential cross sections for the yp —> 7r°p reaction have been measured with the A2 tagged- 
photon facilities at the Mainz Microtron, MAMI C, up to the center-of-mass energy IF = 1.9 GeV. 

The new results, obtained with a fine energy and angular binning, increase the existing quantity of 
7r° photoproduction data by ~ 47%. Owing to the unprecedented statistical accuracy and the full 
angular coverage, the results are sensitive to high partial-wave amplitudes. This is demonstrated by 
the decomposition of the differential cross sections in terms of Legendre polynomials and by further 
comparison to model predictions. A new solution of the SAID partial-wave analysis obtained after 
adding the new data into the fit is presented. 

PACS numbers: 12.38.Aw,13.60.Rj,14.20.-c,25.20.Lj 


I. INTRODUCTION 


‘corresponding author, e-mail: prakhov@ucla.edu Measurements of pion photoproduction on both pro- 

t corresponding author, e-mail: igor@gwu.edu ton and quasifree-neutron targets have a very long his¬ 

tory, started about 65 years ago. Using the first 
bremsstrahlung facilities, pioneering results for yp —> 
7r°p [ll], 7 p —> 7r + n 0, and yn —> tt~p [.T were ob- 
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tained. Despite all the shortcomings of the first measure¬ 
ments (such as large normalization uncertainties, wide 
energy and angular binning, limited angular coverage, 
etc.), those data were crucial for the discovery of the 
first excited nucleon state, A(1232)3/2 + Q. 

Though present electromagnetic facilities in combina¬ 
tion with modern experimental detectors have allowed 
a significant improvement in the quantity and quality 
of pion-photoproduction data Q, full understanding of 
pion-photoproduction dynamics is far from established, 
even in the best investigated area, comprising the first, 
second, and third resonance regions. The properties of 
many nucleon states in this_energy range, especially in 
the third resonance region [fjj, are still not well under¬ 
stood. Resolving the complex partial-wave structure of 
pion photoproduction with minimal model assumptions 
requires more precise data for differential cross sections, 
with fine energy binning and full angular coverage, in 
combination with measurements of various polarization 
observables. 

First results for polarization degrees of freedom have 
been recently reported by the A2, CBELSA/TAPS, and 
GRAAL Collaborations. The beam asymmetry E has 
been measured by GRAAL [t| and CBELSA/TAPS Q. 
Data for the other two single polarization observables T 
(target polarization) and P (recoil-nucleon polarization 
measured as a double-polarization observable with trans¬ 
versely polarized target and linearly polarized beam) 
have been reported from CBELSA/TAPS 0). First re¬ 
sults for double-polarization asymmetries have also been 
measured at CBELSA/TAPS with a longitudinally polar¬ 
ized target (observable G with linearly polarized beam Q 
and observable E with circularly polarized beam [10|) 
and with a transversely polarized target and linearly po¬ 
larized beam (observable H i)- The A2 Collaboration 
at MAMI obtained results for the double-polarization 
asymmetry C* using a circularly polarized photon beam 
and a nucleon recoil polarimeter if]. Further data sets 
are currently under analysis, so that one can expect new 
analyses which provide much tighter constraints on the 
7 p —> 7 r°p reaction in the nearest future. 

This work contributes to these efforts by presenting a 
new high-statistics measurement of the yp —> n°p differ¬ 
ential cross sections conducted by the A2 Collaboration 
for incident-photon energies, 2£ 7 , from 218 MeV up to 
1573 MeV [or center-of-mass (c.m.) energies W =1136- 
1957 MeV]. The data are obtained with a fine binning in 
Ery (~ 4 MeV for all energies below E 1 = 1120 MeV) and 
30 angular bins, covering the full range of the tt° produc¬ 
tion angle. The data obtained above E 1 = 1443 MeV 
{W = 1894 MeV) have limited angular coverage. 

A more detailed analysis of the present yp —> 7 r°p dif¬ 
ferential cross sections, combined with results from the 
A2 Collaboration obtained by measuring polarization ob¬ 
servables, is currently in progress and will be published 
separately. 


II. EXPERIMENTAL SETUP 

The reaction yp —> 7 r°p was measured using the Crystal 
Ball (CB) [lH as a central calorimeter and TAPS fl3l.fl4l| 
as a forward calorimeter. These detectors were installed 
at the energy-tagged bremsstrahlung-photon beam pro¬ 
duced from the electron beam of the Mainz Microtron 
(MAMI) 0,(3. 

The CB detector is a sphere consisting of 672 optically 
isolated Nal(Tl) crystals, shaped as truncated triangular 
pyramids, which point toward the center of the sphere. 
The crystals are arranged in two hemispheres that cover 
93% of 47 t sr, sitting outside a central spherical cavity 
with a radius of 25 cm, which is designed to hold the 
target and inner detectors. 

In the A2 experiments at MAMI C, TAPS was initially 
arranged in a plane consisting of 384 BaF 2 counters of 
hexagonal cross section. It was installed 1.5 nr down¬ 
stream of the CB center and covered the full azimuthal 
range for polar angles from 1° to 20°. In the present ex¬ 
periments, 18 BaF 2 crystals, covering polar angles from 
1° to 5°, were replaced with 72 PbW 04 crystals. This 
allowed running with a much higher MAMI electron cur¬ 
rent, without decreasing the TAPS efficiency due to the 
very high count rate in the crystals near the photon- 
beam line. The energy resolution of the PbW 04 crystals 
has not been understood well yet, and the use of their 
information was restricted in the present analysis, as de¬ 
scribed later. More details on the calorimeters and their 
resolutions are given in Ref. E3 and references therein. 

The present measurements used 1557-MeV and 1604- 
MeV electron beams from the Mainz Microtron, MAMI 
C |l 6 ]. The data with the 1557-MeV beam were taken 
in the first half of 2013 (Run-I) and those with the 
1604-MeV beam in the second half of 2014 (Run-II). 
Bremsstrahlung photons, produced by the beam elec¬ 
trons in a lO-pm Cu radiator and collimated by a 4- 
mm-diameter Pb collimator, were incident on a 10-cm- 
long liquid hydrogen (LH 2 ) target located in the center 
of the CB. In Run-I with the 1557-MeV electron beam, 
the energies of the incident photons were analyzed up 
to 1448 MeV by detecting the postbremsstrahlung elec¬ 
trons in the Glasgow tagged-photon spectrometer (Glas¬ 
gow tagger) Il8l - [20f . The uncertainty in the energy of 
the tagged photons is mainly determined by the number 
of tagger focal-plane detectors in combination with the 
energy of the MAMI electron beam used in the experi¬ 
ments. Increasing the MAMI energy results in increasing 
both the energy range covered by the spectrometer and 
the corresponding uncertainty in E 1 . For the MAMI en¬ 
ergy of 1557 MeV and the Glasgow tagger, such an un¬ 
certainty was about ±2 MeV. The systematic uncertainty 
in the absolute value of £/,, which is dominated by the 
energy calibration of the tagger, was about 0.5 MeV [20| . 
More details on the tagger energy calibration and uncer¬ 
tainties in the energies can be found in Ref. [2l[ . 

In Run-II with the 1604-MeV electron beam, the 
energies of the incident photons were analyzed from 
1426 MeV up to 1576 MeV by detecting the post- 
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FIG. 1: cos 9 distributions for 7 p —¥ n°p, with 9 being the angle between the directions of the outgoing n° and the incident 
photon in the c.m. frame, obtained for E 7 = 318 MeV (W = 1215 MeV): (a) experimental data after subtraction of the 
empty-target and the random backgrounds; (b) angular acceptance obtained from the MC simulation of 7 p —> n °p; (c) results 
for the 7 p —>• n°p differential cross section compared to the prediction by SAID [2J| with its CM12 [24] solution. The error bars 
on all data points represent statistical uncertainties only. 



bremsstrahlung electrons in the Mainz end-point tagger 
(EPT) with 47 focal-plane detectors. Since the Glasgow 
tagger was not designed to measure the high-energy tail 
of the bremsstrahlung spectrum, the EPT spectrometer 
was built to conduct 7 / measurements by covering this 
low-energy range of postbremsstrahlung electrons. The 
uncertainty of the EPT in E 1 due to the width of its 
focal-plane detectors was about ±1.6 MeV, with a simi¬ 
lar value (i.e., ~ 1.6 MeV) in the systematic uncertainty 
in Ery due to the EPT energy calibration. The energy 
calibration of the EPT is based only on the simulation 
of electron tracing, using measured magnetic-field maps. 
The correctness of this calibration, as well as its uncer¬ 
tainty, was checked by measuring the position of the rf 
threshold, E 1 ~ 1447 MeV. 

Because the main goal of Run-I was to measure the 
7 p — > 7 r°p reaction from E.y = 218 MeV up to 1448 MeV 
and with full coverage of the ir° production angle, the 
trigger required the total energy deposited in the CB 
to exceed ~120 MeV. Run-II was mainly dedicated to 
studying if physics. For that, the experimental trigger 
required the CB energy to exceed ~540 MeV, rejecting 
the detection of 7r° production in the very forward and 
backward angles. 


III. DATA HANDLING 

The reaction qp — > ir°p — > qqp was searched for in 
events reconstructed with two or three clusters detected 
in the CB and TAPS together. In the analysis of the data 
from Run-I, the information from the 72 PbW 04 crys¬ 
tals of TAPS was not used at all. In the analysis of the 
data from Run-II, in which all PbW 04 crystals demon¬ 
strated good performance, their information was used in 
the cluster reconstruction. The two-cluster events were 
analyzed assuming that the final-state proton was not 
detected. This typically happens when the outgoing pro¬ 
ton is stopped in the material of the downstream beam 
tunnel of the CB, or the proton kinetic energy in the lab¬ 
oratory system is below the software cluster threshold, 
which was 12 MeV in the present analysis. 

The selection of event candidates and the recon¬ 
struction of the reaction kinematics were based on the 
kinematic-fit technique. Details on the kinematic-fit 
parametrization of the detector information and resolu¬ 
tions are given in Ref. H3- The kinematic-fit technique 
was also used for the offline calibration of the calorime¬ 
ters. The resolution functions used in the kinematic fit 
were adjusted to result in proper stretch functions (or 








4 



FIG. 3: Same as Fig. [T] but for E 7 = 1512 MeV (W = 1928 MeV). 


pulls) and probability distributions. All two- and three- 
cluster events that satisfied the yp —> n°p — > yyp hypoth¬ 
esis with a probability greater than 1 % were accepted 
for further analysis. The kinematic-fit output was used 
to reconstruct the four-momenta of the outgoing parti¬ 
cles. The events from the data of Run-II with more than 
one cluster in the 72 PbWC >4 crystals of TAPS were dis¬ 
carded from the analysis, and the events with one clus¬ 
ter in these crystals were tested only for the hypothesis 
assuming that this cluster is produced by the recoil pro¬ 
ton. Such a restriction was imposed due to the unknown 
energy-resolution function for the PbWC >4 crystals, and 
only the angular information from the proton cluster was 
used in the kinematic fit. 

The determination of the experimental acceptance was 
based on a Monte Carlo (MC) simulation of yp —> 7r°p —> 
77 p with an isotropic production-angular distribution 
and a uniform beam distribution generated for the full 
energy range of the tagged bremsstrahlung photons. All 
MC events were propagated through a GEANT (3.21) 
simulation of the experimental setup. To reproduce the 
experimental resolutions, the GEANT output was sub¬ 
ject to additional smearing, thus allowing both the simu¬ 
lated and experimental data to be analyzed in the same 
way. This additional smearing was adjusted by describ¬ 
ing the experimental invariant-mass resolutions as well as 
the kinematic-fit probability distributions. The final ad¬ 
justment of the detection efficiency for yp —>• 7r°p —> yyp 
took into account the trigger requirements in the analysis 
of the MC events. 

MC simulations for other possible reactions (such as 
yp —> rjp — > yyp, yp — >• 7r°7r°p, yp — >• uip — >• 7r°yp, 

yp —> ir°ir + n) showed that background contributions 
from these sources were negligibly small. Thus, exper¬ 
imental spectra with selected events were contaminated 
only by two sources of background events, the distribu¬ 
tions of which were directly subtracted afterwards. The 
first source of background was due to interactions of in¬ 
cident photons in the windows of the target cell. The 
subtraction of this background was based on the analy¬ 
sis of data samples that were taken with an empty (no 
liquid hydrogen) target. The weight for the subtraction 
of the empty-target spectra was taken as a ratio of the 


photon-beam fluxes for the data samples with the full 
and the empty target. A second background was caused 
by random coincidences of the tagger counts with the 
experimental trigger; its subtraction was carried out by 
using event samples for which all coincidences were ran¬ 
dom (see Refs. [ 13,1221 for more details). 

For measuring the yp —> 7r°p differential cross sections, 
all selected events were divided into 30 equal-width cos 9 
bins, covering the full range from -1 to 1, where 9 is the 
angle between the directions of the outgoing 7r° and the 
incident photon in the c.m. frame. Also, the quantity of 
data collected for yp —>• 7r°p was sufficient to obtain sta¬ 
tistically accurate differential cross sections for every tag¬ 
ger channel in the energy range from F 7 = 218 MeV to 
1120 MeV, providing ~ 4 MeV binning in J5 7 . For higher 
energies, energy bins were combined in two or more tag¬ 
ger channels. 

The typical experimental statistics and acceptance are 
illustrated in Figs. Ill [2j and [3] for different energies and 
data sets. Figure Oja) shows the experimental cos 9 dis¬ 
tribution for F 7 = 318 MeV (W = 1215 MeV), which 
is obtained after subtracting the random and the empty- 
target backgrounds. Since this energy bin is close to the 
A(1232)3/2 + maximum, the statistical uncertainties in 
this distribution, based on 8 x 10 6 events, are very small 
(~ 0.2%). Figure [D(b) shows the corresponding angu¬ 
lar acceptance, which is greater than 80% for the central 
range of 7r° production angles. The resulting yp —> 7r°p 
differential cross section for E 1 = 318 MeV is depicted 
in Fig. [ljc). It is also compared to the prediction from 
the partial-wave analysis (PWA) by SAID [23[ with its 
CM12 [ 24 } solution (shown by a solid line). The experi¬ 
mental data points lie slightly above the phenomenologi¬ 
cal prediction, but their angular dependence is very close 
to it. Because all modern PWAs give very similar pre¬ 
dictions for the yp —> 1 r°p differential cross section at the 
A(1232)3/2 + maximum, the present results in this region 
do not contradict any of them. 

Figure HJa) shows the experimental cos 9 distribution 
for E 1 = 1202 MeV {W = 1771 MeV), which is near the 
third resonance region. The statistical uncertainties in 
this distribution, based on 8.5 x 10 4 events, are larger but 
are still sufficient for a reliable analysis of the resulting 
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FIG. 4: (Color online) Comparison of the 7 p —> 7r°p differential cross sections obtained by the A2 Collaboration in different 
experiments: (a) results for E 7 = 218 MeV from Run-I (black points) and from the experiment [25| with MAMI energy of 
855 MeV (green points); (b) same as (a) but for E~, = 331 MeV; (c) results for E 7 = 1434 MeV from Run-I (black points) 
and Run-II (red points); (d) same as (c) but for E 7 = 1443 and 1441 MeV of Run-I and Run-II, respectively; The error bars 
on all data points represent statistical uncertainties only. 


differential cross section. Since the shape of the 7 p —> n°p 
differential cross sections in this region changes rapidly 
with energy, the main goal of the present work was to 
measure it with the finest energy binning. Figure [2fb) 
shows that the angular acceptance drops for the very 
forward angles. At these energies, this occurs as many it 0 
mesons produced at forward angles have both their decay 
photons going into TAPS, and such events do not pass 
the trigger requirements. The resulting differential cross 
section, plotted in Fig. 0(c), demonstrates again a quite 
reasonable agreement with the SAID CM12 solution 24], 
Figure [3] shows the cos 9 distributions for E 1 = 
1512 MeV (W = 1928 MeV) from the data taken in Run- 
II. Despite the large quantity of the experimental events 
(9 x 10 5 ), a very low, and even zero, acceptance for the 
very forward and backward angles, caused by the trigger, 
did not allow the first two and the last two cos 9 bins to 
be measured. However, the remaining 26 data points of 
the differential cross section lie close to the prediction 
from SAID CM12 [H|. 


IV. EXPERIMENTAL RESULTS AND THEIR 
DISCUSSION 

The yp —1 7 r°p differential cross sections were obtained 
by taking into account the number of protons in the LH 2 
target and the photon-beam flux from the tagging fa¬ 
cilities corrected for the fraction rejected by collimators. 
The overall systematic uncertainty due to the calculation 
of the detection efficiency and the photon-beam flux was 
estimated to be 4% for the data taken in Run-I and 5% 
for the data taken in Run-II. 

Another source of systematic uncertainty is due to im¬ 
perfections in reproducing the angular dependence of the 
cos 9 acceptance. Such imperfections are easier to see 
in differential cross sections with very low statistical un¬ 
certainties and a smooth shape. In the present data, it 
can be seen, for example, on the top of the experimental 


spectrum shown in Fig. [He). The angular dependence of 
differential cross sections can also be smeared due to the 
limited resolution in the n° production angle, resulting 
in a larger population of angular bins with a lower cross 
section. Such an effect is, possibly, seen in Fig. (He) near 
cos 9 = 0.5. 

Possible magnitudes of systematic uncertainties in dif¬ 
ferential cross sections can be checked by comparing sim¬ 
ilar results obtained by analyzing different data. Figure[2 
illustrates self-consistency of the results for the 7 p —> n°p 
differential cross sections obtained by the A2 Collabora¬ 
tion in three different experiments. The results for the 
lowest energy, E 1 = 218 MeV, taken in Run-I are com¬ 
pared in Fig. [2a) with similar results obtained in the 
A2 experiment (A2-2008) with the Glasgow tagger and 
MAMI energy of 855 MeV (25|. As seen, the agreement 
of the two independent measurements is fairly good for 
both the absolute values and the angular dependence. 
Figure [4jb) compares the results of Run-I and A2-2008 
for an energy near the A(1232)3/2" t ' maximum, showing a 
small discrepancy in the angular dependence. Continuing 
the comparison of the results from Run-I and A2-2008 at 
energies above the A(1232)3/2 + peak makes little sense 
as, at those energies, the A2-2008 data have much larger 
systematic uncertainties in the photon flux and in the 
angular acceptance, affected by the trigger. The results 
for two the highest energies, E 1 = 1434 and 1443 MeV, 
taken in Run-I are compared in Figs. [2c) and (d) with 
similar results from the present work obtained in Run-II. 
Despite a general agreement, there are some discrepan¬ 
cies between the two sets of results, which cannot be 
explained by their statistical uncertainties. In our opin¬ 
ion, such discrepancies could be caused by several factors. 
There was a significant difference in the trigger of the two 
experiments and in the energy resolution of the calorime¬ 
ters, as the experiment with the EPT (Run-II) was made 
with a much higher MAMI electron current. A poorer en¬ 
ergy resolution in Run-II resulted in a larger smearing of 
the cos 9 angular dependence. Another possible explana- 
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FIG. 5: (Color online) Comparison of the -yp —> n°p differential cross sections from the present work (black points) with 
most recent results from other experimental setups at MAMI (blue points) and other facilities: GRAAL 0] (red points), 
CL AS [27 ] (green points), CB-ELSA [2Sj (magenta points), and CBELS A/TAPS [29| (cyan points). Values of E-, in each panel 
indicate the photon-beam energies for the results of this work. All data from other experiments are shown for nearby energies. 
The error bars on the data points from this work represent statistical uncertainties only. The uncertainties of some previous 
results also include angular-dependent systematic uncertainties. 


tion is a slight mismatch in the energy calibrations of the 
two different tagged-photon spectrometers. Overall, the 
average systematic uncertainty in the differential cross 
sections due to the resolution in the c.m. 9 and possible 
imperfections in the angular acceptance was estimated to 
be 2% for the results of Run-I. The corresponding sys¬ 
tematic uncertainty for the results of Run-II is at least 
twice as high. 

Comparison of the results for the yp —> 7r°p differential 
cross sections from the present work with the most recent 
published data from other collaborations l2bl429[ is il¬ 
lustrated in Fig. [5] for eight different energies. Previous 
results with either large angular bins or very restricted 
angular coverage are not included in the comparison. As 
seen in Fig. [5j data with full angular coverage existed 
only in the region below the third resonance region. The 
results of this work are in a general agreement with all the 
data shown, although those data are not always entirely 
consistent with each other. As also evident, the present 
data from Run-I have much superior statistical accuracy, 
combined with a finer angular and energy binning. 

The total cross sections for yp —> n°p were determined 
by integrating the differential cross sections for the data 
obtained from Run-I, which have the full angular cov¬ 
erage. The resulting total cross sections are plotted in 
Fig. [6] as a function of the c.m. energy W along with ear¬ 
lier A2 results obtained for energies up to = 444 MeV 
(W = 1309 MeV) [25| and results from Ref. [26j, obtained 
for energies up to E 1 = 790 MeV {W = 1537 MeV), and 


from Ref. [28} , obtained for energies above E 7 = 300 MeV 
(W = 1201 MeV). The latter data have no full angular 
coverage; they were obtained with extrapolating experi¬ 
mental differential cross sections into forward and back¬ 
ward angles by using the corresponding PWA results. For 
completeness, the data at higher energies are also com¬ 
pared in Fig. EDb) to the SAID solution CM12 [24|], shown 
by a magenta dashed line. As seen in Fig. [^a) and (b), 
the present results are in good agreement, within the to¬ 
tal uncertainties (calculated by adding the statistical and 
the overall systematic uncertainty in quadrature), with 
the earlier A2 results and in general agreement with the 
most recent experimental data and SAID CM2 solution 
based on htting previous data, even if they do not have 
the full angular coverage. 


Comparison of the present results for the yp —> ir°p 
differential cross sections with predictions of different 
partial-wave and coupled-channel analyses is illustrated 
for selected energies in Fig. [7] Three sets of predic¬ 
tions shown in this figure correspond to PWAs (namely, 
SAID CM12 Jp, MAID2007 0, and Bonn-Gatchina 
BG2014-02 |3lj) based on fitting previous data. The 
fourth set of curves comes from SAID PR15 solution 
obtained after adding the present data points into the 
SAID fit. Compared to the SAID solution CM12, the 
new fit PR15 also involves the data added to obtain SAID 
DU13 solution [ 5 } and other most recent polarization data 
from CBELSA/TAPS @,U3l and A2 at MAMI [Tl(. As 
seen, all three the previous PWAs are very close to each 
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FIG. 6 : (Color online) Total cross sections for 7 p — ► tt°p from this work (black points) as a function of the c.m. energy W 
compared with earlier A2 results [25j (red points), with results from Refs. [2f| (blue points) and [28|] (green points), and with 
the SAID solution CM12 [24lj shown in (b) by a magenta dashed line. For a better comparison, the energy range is divided 
in two intervals: (a) W < 1340 MeV and (b) W > 1340 MeV. The error bars on the experimental data points represent their 
total uncertainties. 


TABLE I: Average x 2 /dp values (including information on the total y 2 from Eq. JTJ and the number of data points used) for 
MAID2007 [3C| (fitting data only below W = 2 GeV) and three SAID solutions: CM12 [24|, DU13 and PR15. 


X 2 /dp 

MAID2007 

SAID CM12 

SAID DU13 

SAID PR15 

present yp —> n°p 

99722/7978 = 12.5 

24036/7978 = 3.0 

28745/7978 = 3.6 

9434/7978 = 1.2 

previous 

7 p —> 7 r°p 

7 p 7 r + n 

7 n —► 7 r~p 

7 n —> 7 r°n 

135195/15468 = 8.7 
171853/8092 = 21.2 
25335/2806 = 9.0 
22568/364 = 62.0 

73211/17087 = 4.3 
23533/8959 = 2.6 
53657/3162 = 17.0 
979/364 = 2.7 

41500/17087 = 2.4 
19312/8959 = 2.2 
7561/3162 = 2.4 
1436/364 = 4.0 

40793/17087 = 2.4 
17540/8959 = 2.0 
6572/3162 = 2.1 
1187/364 = 3.3 


other and to the present data points in the energy re¬ 
gion around the A(1232)3/2 + maximum. Differences 
between the PWAs themselves as well as between the 
PWAs and the data points increase, especially for the 
MAID2007 solution, at energies that are away from the 
A(1232)3/2 + peak. Then all the PWAs and the data 
points become closer again in the second and the third 
resonance regions. The largest difference between the 
present data and the three previous PWAs is observed 
for the MAID2007 solution; this could be explained by 
the fact that many of the most recent data sets were not 
available at the time when the MAID2007 solution was 
released. 

Quantitatively, the agreement of the present data with 
PWAs can usually be estimated via y 2 values calculated 
from deviations of experimental data points from PWA 
solutions. The y 2 function used in the minimization pro¬ 
cedure by SAID is given by formula Q 


X 


2 


y^X9.i - 6| XP ) 2 + - 1 


(1) 


where 0 ° xp is an individual experimental value of the 


differential cross section with its uncertainty ti (includes 
both the statistical and the individual systematic uncer¬ 
tainty discussed in the text above), 0 ; is the fit value 
calculated for the same energy and the production an¬ 
gle, and A is a normalization parameter with its uncer¬ 
tainty ex that is determined by the overall systematic 
uncertainties of the data (see start of this section). The 
average y 2 /dp values, along with the information on the 
total y 2 from Eq. 0 and the number of da ta p oints used 
(dp), are listed in Table Q] for MAID2007 [30j and three 
SAID solutions: CM12 Hi, DU13 @, and PR15. As 
seen, the numbers in Table U confirm that the average 
y 2 /dp is worse for the older solution MAID2007 than for 
the more recent SAID solutions CM12 and DU13. The 
magnitudes of y 2 /dp obtained for the new SAID solution 
PR15 indicate that the new yp —>• 7 r°p data are consis¬ 
tent with the existing 7 N —> ttN data sets. The number 
of data points available for each photoproduction reac¬ 
tion, which are listed in Table HI shows that the present 
7 p —> 7 r°p data increase the existing yp —> 7r°p statis¬ 
tics by ~ 47% and the existing yiV —>• 7 tN statistics by 
~ 27%. 

The y 2 /dp values, averaged within each energy bin, 


















FIG. 7: (Color online) Selected results of this work (black points) for the 7 p —> n°p differential cross sections compared to 
existing PWA solutions from SAID CM12 [24[ (blue dashed line), MAID2007 13C3J (magenta long-dashed line), and Bonn- 
Gatchina BG2014-02 [3l|] (green dash-dotted line) and to a new SAID PR15 solution (red solid line) obtained after adding the 
present data points into the fit. The error bars on all data points represent statistical uncertainties only. 


are plotted as a function of E 1 in Fig. Ufa) for the previ¬ 
ous SAID solution DU13 [H and for the new SAID PR15 
solution, obtained after adding the present data into the 
fit. This comparison demonstrates the energies at which 
the new data give the most significant impact. The re¬ 
sults obtained for the normalization parameter X from 


the DU13 and PR15 solutions are plotted as a function 
of E-y in Fig. (HKb). 

Another useful presentation of the new 7 p —>- n°p re¬ 
sults is obtained by plotting them as an excitation func¬ 
tion at a particular production angle 6 of 7 r°. In Fig. 0 
the differential cross sections are shown as a function of 
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FIG. 8: (Color online) Comparison of the 
previous SAID solution DU13 [|| (blue open 
triangles) applied to the present data and 
the new SAID PR15 solution (red full cir¬ 
cles) obtained after adding the present data 
into the fit for: (a) x 2 /dp values averaged 
only within each energy bin E 1: where the 
horizontal lines (blue dashed for DU13 and 
red solid for PR15) show the corresponding 
X 2 /dp values from Table U (b) results for 
the normalization parameter X as a func¬ 
tion of E-y, where the horizontal solid lines 
show the deviations (from the dashed line 
at X = 1) equal to the magnitude of the 
overall systematic uncertainties, which are 
4% and 5% for the data from Run-I and 
Run-II, respectively. 


the c.m. energy W for ten of 30 measured cos 9 bins. As 
seen, the peaks corresponding to the first, the second, and 
the third resonance region are clearly evident in the exci¬ 
tation function for each n° production angle. The exper¬ 
imental distributions are also compared to the three pre¬ 
vious PWA solutions (SAID CM12 [3, MAID2007 l30l| - 
and BG2014-02 [§T|) and the new SAID PR15 solution. 
For a better comparison in the energy range away from 
the A(1232)3/2 + peak, all distributions are also plotted 
after rescaling by a factor of 5. The comparison demon¬ 
strates that the largest discrepancies are seen for the very 
forward and backward production angles. 

The strongest dip in the yp —> ir°p excitation func¬ 
tion is observed at very backward 7r° angles at the c.m. 
energies slightly below W = 1.4 GeV. An example of 
the yp —> 7r°p differential cross section near this dip 
(namely, at W = 1380 MeV) is shown in Fig. flOl'a). 
This effect, caused by the nearly complete cancellation 
between different contributions to the production cross 
section, is even better seen in the results of the new SAID 
PR15 solution, shown in Fig. fTtilb) by a blue solid line 
for cos# = —1 (or 9 = 180°). It is interesting that a 
very similar cancellation in the very backward 7r° produc¬ 
tion angles was observed in the charge-exchange reaction 
7 T~p —> 7 r°n at the same energy range. The SAID WI08 
solution |33] (rescaled by a factor 10 -3 ) for 7r~p —> n°n 
at 9 = 180° is plotted in Fig. HHIb) by a red dashed 
line. In contrast, the yp —j tt + n reaction, depicted for 
the SAID PR15 solution by a black dashed line, has no 
such feature. Unfortunately, the quality and the angular 
coverage of the existing y n —> tt° n data do not allow a re¬ 
liable prediction for very backward production 7r° angles 
in this energy range [34| . 

For the c.m. energies above W = 1.4 GeV, the present 
results for very backward 7r° angles are in good agreement 
with similar results ( 9 = 170°) published in Ref. 3rl . 
also showing a cusp structure at the 77 -photoproduction 
threshold {W = 1.486 GeV). 

The statistical accuracy of the present results is suffi¬ 
cient for reliable fits of the yp —> 7r°p differential cross 
sections with series of Legendre polynomials up to a high 


order. Of course, such decompositions are much less re¬ 
liable without full angular coverage of the 7r° produc¬ 
tion angle. Thus, the present results from Run-II should 
be used only in energy-dependent PWAs with additional 
model assumptions. For Run-I, we describe the yp —> 7 r°p 
differential cross sections of each energy bin in terms of 
a series of orthogonal Legendre polynomials Pj (cos 9): 

7 Jmax 

^(W, cos 9) = Y M w ) Pj {cos 9), 
j =0 

with the energy-dependent decomposition coefficients 

MW). 

From the phenomenological point of view, the number 
of terms required in the series, j m ax> is related to the 
angular momenta of essential partial-wave amplitudes or 
to the highest spin of resonances, existing in the studied 
energy range. An isolated resonance with spin J = [2i + 
l)/2 can contribute to coefficients Aj only with even j 
up to j = 2 £ and, hypothetically, to other coefficients 
via interference with sizable background amplitudes. For 
example, a wide four-star resonance A (2420) with J p = 
11/2+ [6], although lying above the energy range of the 
present measurements, can affect Aj up to j = 10 by 
itself and even higher Aj via the interference with other 
states. On the other hand, the maximum order j ma x, up 
to which such a decomposition is meaningful, is limited 
by the data quality. 

In this work, the data of Run-I were fitted up to or¬ 
der jmax = 10. It was observed that including Legendre 
polynomials beyond j m ax = 8 in the fits could not sig¬ 
nificantly improve the reduced x' 2 (he., divided by the 
number of degrees of freedom, \' 2 /ndf). The comparison 
of the reduced x 2 for the fits of the data from Run-I with 
jmax = 6, 8, and 10 is shown in Fig. [TTJ using only the 
statistical uncertainties. In the energy range around the 
A(1232)3/2 + region, which is below W = 1.3 GeV, fit¬ 
ting with jmax = 6 should, in principle, be sufficient for 
a fairly good description of the differential cross sections. 
However, the high statistical accuracy, which arises from 
the very large cross sections, results in the reduced \ 2 
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FIG. 9: (Color online) Present results for the 7 p —» n°p excitation function as a function of the c.m. energy W, shown for ten 
of 30 measured cos 8 bins (where 8 is the 7 r° production angle), are compared to previous PWA solutions from SAID CM12 [24]] 
(blue dashed line), MAID2007 [3C| (magenta long-dashed line), and BG2014-02 [3lJ (green dash-dotted line) and to a new SAID 
PR15 solution (red solid line). The error bars on all data points (black points for Run-I and cyan points for Run-II) represent 
statistical uncertainties only. For a better comparison, all distributions are also plotted after rescaling by a factor of 5. 


values varying mostly between two and three. As also 
seen in Fig. 1111 including higher-order coefficients in the 
fits does not improve the situation much. It is necessary 
to take into account angular-dependent systematic uncer¬ 
tainties, with values varying between 1% and 2%, in order 


to obtain reasonable y 2 values in the A(1232)3/2 + re¬ 
gion, even in the fits with j max = 6. Above W =1.7 GeV, 
using j max = 6 is not sufficient, and the data require the 
minimum order of j max = 8, corresponding to G -wave 
amplitudes. 
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FIG. 10: (Color online) (a) present 
results for the 7 p — > -K°p differential 
cross section (blue full circles) at E 7 = 
546 MeV (W = 1380 MeV) fitted with 
Legendre polynomials (red dashed line); 
(b) the new SAID PR15 solution (blue 
solid line) for n° production angle 9 = 
180° (cos# = — 1 ) as a function of the 
c.m. energy W, compared with a simi¬ 
lar n~p —> 7 r°n prediction (rescaled by a 
factor 10“ 3 ) from SAID WI08 [H (red 
dashed line) and a -yp — > -K + n predic¬ 
tion (rescaled by a factor 0.3) from SAID 
PR15 (black dashed line). 



FIG. 11: (Color online) Energy dependence of the reduced 
X 2 values for Legendre fits of the present 7 p —> n°p differ¬ 
ential cross sections with j max = 6 (black triangles down), 8 
(blue circles), and 10 (red triangles up), where the angular- 
dependent systematic uncertainties were not taken into ac¬ 
count. 


The results of the Legendre-polynomial fits for each 
coefficient Aj are depicted in Fig. [T2l showing their en¬ 
ergy dependence in unprecedented resolution. The Aj 
error bars represent uncertainties obtained from fits that 
took into account only the statistical uncertainties in the 
cross sections. As previously discussed, the effect from 
the angular-dependent systematic uncertainties should 
be considered, especially for the A(1232)3/2 + energy re¬ 
gion. To estimate the impact from these uncertainties, 
the values of the differential cross section for each an¬ 
gular bin were randomly shifted within the uncertainty 
magnitude (2%), followed by a new Legendre-polynomial 
fit. Such a procedure was then repeated several times, 
and the average spread of the fit results for each Aj was 
considered as its systematic uncertainty caused by the 
angular-dependent systematic uncertainties in the differ¬ 
ential cross sections. The systematic uncertainties in Aj, 
obtained from the fits with j max = 10 , are depicted by 
a black histogram at the bottom of each plot in Fig. [12J 
clearly peaking at the position of A(1232)3/2 + . In this 
region, the magnitudes of the systematic uncertainties 
in Aj at high orders are very close to the magnitudes 
of Aj themselves. The Legendre-polynomial fits for the 


c.m. energies below W = 1280 MeV were made for both 
jmax = 6 and j ma x = 10. Their comparison showed good 
agreement of the results for Aq to Aq. 

One of the most obvious features in the energy de¬ 
pendence of Aj is the absence of any well-defined sharp 
structures, which could be associated with unknown nar¬ 
row resonances. The results for A 0 show good agree¬ 
ment with the values for the total cross sections obtained 
by the integration of the differential cross sections, pro¬ 
viding another confirmation of the good quality of the 
present data. The Aq energy dependence clearly demon¬ 
strates the first, second, and third resonance regions, 
but appears structureless at higher energies. In con¬ 
trast, some of the other Aj reveal structures even at 
the energies above the third resonance region. These 
structures could be associated with contributions from 
N* and/or A* states, with spins up to 7/2 or higher, 
that are known in this region or near it 6 |. These 
resonances have masses close to the upper end of the 
studied energy range, as A(1950)7/2 + , or even above it, 
as 77(2190)7/2“ N( 2220)9/2+, and 77(2250)9/2“, but 
their large widths (~ 300 MeV or even larger) make those 
resonance contributions noticeable in the present data. 

In Fig. [T2j the coefficients Aj are also compared to re¬ 
sults obtained from the predictions generated with SAID 
CM12 [H, MAID2007 fl|, BG2014-02 [31.], and SAID 
PR15, by fitting them, after adding small uncertainties, 
with Legendre polynomials up to j max = 10. As seen, the 
Legendre-polynomial fits to the predictions of all mod¬ 
els reproduce fairly well only the coefficients Aq and A 2 . 
Partial agreement is observed for A\ and the coefficients 
A 3 to A 7 . The coefficients Ag to A 10 are small and show 
some meaningful structures only above W = 1.7 GeV. 

Of course, such a Legendre decomposition of the 
present differential cross sections is only a very first step 
towards a full single-energy partial-wave analysis. A di¬ 
rect extraction of partial-wave amplitudes by simulta¬ 
neous fitting the present differential cross sections and 
polarization observables will be performed in the future. 
However, the present analysis already shows the sensitiv¬ 
ity of the new data to partial-wave amplitudes up to G 
waves (£ = 4). 
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FIG. 12: (Color online) Coefficients Aj (in pb/sr units) obtained from fitting the present 7 p —> n°p differential cross sections 
with Legendre polynomials up to order six (cyan points) and ten (green points), compared to the results of the Legendre- 
polynomial fits ()max = 10) made to the predictions generated from SAID CM12 [juj (blue solid line), MAID2007 [30] (red 
solid line), BG2014-02 [HI] (black solid line), and SAID PR15 (blue dashed line). The error bars on all points represent Aj 
uncertainties from the fits with using only the statistical uncertainties. A black histogram at the bottom of each plot shows 
the systematic uncertainty in Aj caused by the angular-dependent systematics in the differential cross sections. 


V. SUMMARY AND CONCLUSIONS 

New measurements of n° photoproduction on the pro¬ 
ton have been conducted with the A2 tagged-photon 
facilities at energies provided by the Mainz Microtron, 
MAMI C, with a maximum electron-beam energy of 
1.6 GeV. The differential cross sections, obtained with 
a fine energy and angular binning, increase the existing 
quantity of 7r° photoproduction data by ~ 47%. A new 
solution PR15 of the SAID partial-wave analysis obtained 
after adding the new data into the fit is presented. The 
magnitudes of y 2 /dp obtained for this solution indicate 
that the new jp —»• n°p data are consistent with the ex¬ 
isting 7 N —> 7 tN data sets. Owing to the unprecedented 
statistical accuracy and the full angular coverage, the re¬ 
sults are sensitive to high partial-wave amplitudes. This 
is demonstrated by the decomposition of the differential 
cross sections in terms of Legendre polynomials and by 
further comparison to model predictions. The present 
data are expected to be invaluable for future partial-wave 
and coupled-channel analyses, which could provide much 
stronger constraints on the properties of nucleon states 
known in this energy range and, perhaps, even reveal 


new resonances. A more detailed analysis of the present 
7 P —> ir°p differential cross sections, combined with re¬ 
sults from the A2 Collaboration obtained by measuring 
polarization observables, is currently in progress and will 
be published separately. 
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